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The experimentally determined size-evolutionary pat- 
terns (SEPs) of simple metal clusters pertaining to ion- 
ization potentials (IPs), electron afRnities, and monomer 
separation energies have attracted much attention in the 
past ten yearsa'EI, since they may reflect the electronic 
shell structure of clusters (the major features of SEPs 
are associated with steps at magic numbers, but often 
there is in addition substantial fine structure, such as 
odd-even alternations). 

Due to the required computational effort, first- 
principles (FP) theoretical studies (which incorporate the 
ionic geometry) of such SEPs are usually limited to clus- 
ter sizes of the order of ten atomsJj'cl Recently, however, 
systematic theoretical investigations of such SEPs have 
been performed for a broad range of cluster sizes us- 
ing the jallium-related Shell-Correction-Method (SCM) 
approach. Q^B -. 

In a recent article, Gardet et ala have studied the 
IPs of small lithium clusters Li^v (with N < 12) 
using a FP - density functional theory (DFT) and 
found reasonable agreement between theory and exper- 
iment. Furthermore, through a comparison of their re- 
sults with those obtained from Kohn-Sham local-density- 
approximation (K S- | Lp A) calculations on a spherical jel- 
lium background,Eiltll they concluded that the jellium 
model is inadequate for a proper description of the IPs 
of such systems. 

The purpose of this comment is to clarify that, while 
modelling the ions in a cluster via a uniform jellium back- 
ground is certainly an approximation, the above conclu- 
sion of Gardet et al. pertains to limitations introduced 
through neglect of deviations from spherical symmetry, 
rather than to the jellium approximation itself. Indeed, 
in a series of papers,El~El we have demonstrated that con- 
sideration of triaxial shape deformations drastically im- 
proves the agreement between the jellium approximation 
and experiment for all instances of the aforementioned 
SEPs and for sizes up to 100 atoms, as well as for a vari- 
ety of metal species (namely, alkali metals, such and Na 
and K, and noble metals, such as Cu and Ag). 

To further elucidate the importance of shape deforma- 
tions, we display in Fig. 1 the IPs of small Lijv iclusters (in 
the same size-range as with Ref. |9|) calculatedll3 with our 
SCM for three different families of shapes of the jellium 
background [namely, spherical, spheroidal (axially sym- 
metric) , and ellipsoidal (triaxial)]- and compare them to 
the experimental measurements. Ej 

Fig. 1 reveals that spherical shapes (top panel) ex- 
hibit a characteristic sawtoothed profile, well known from 



previous spherical KS-LDA studiest3 and similar to the 
curv(£J labeled jeUium-LDA in Fig. 12 of Ref. ||. Apart 
from major-shell closures, this sawtoothed profile de- 
scribes the data rather poorly (notice in particular the 
absence of fine structure between major shell closures at 
iV = 2, 8 and 20). 

The spheroidal model (middle panel) exhibits substan- 
tial improvement in describing the experimental trend. 
Furthermore, the ellipsoidal case (bottom panel) im- 
proves the agreement between the SCM results and ex- 
periment even further, in particular in the size range 
11 < A^ < 14. The essential improvement introduced 
by the deformations over the spherical case concerns the 
very good description of the subshell closure at A^ = 14 
and of odd-even alternations between major shell clo- 
sures. 

In summary, we have illustrated once again that, 
within the jellium approximation, deformed cluster 
shapes provide an adequate description of the observed 
systematic size dependence of the properties of simple 
metal clusters and should necessarily be employed in 
comparisons with other theoretical approaches. 
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FIG. 1. IPs for Liiv clusters in the range 2 < iV < 16. Solid 
dots: Theoretical results derived from the SCM (see Refs. 
bUtI) using three different shape models for the jellium back- 
ground. Top panel: The spherical model compared to exper- 
imental data (open squares). Middle panel: The spheroidal 
model compared to experimental data. Lower panel: The 
ellipsoidal model compared to experimental data. The exper- 
imental measurements (open squares) were taken from Ref. 
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